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We derive rigorous quantum mechanical bounds for the heat current through a nanojunction 
connecting two thermal baths at different temperatures. Based on exact sum rules, these bounds 
compliment the well-known quantum of thermal conductance kq = -Kk^T/Qh, which provides a 
bound for low-temperature heat transport in all systems, but is saturated only for noninteracting 
transport. In contrast, our bounds are saturated at high temperatures—but still in the quantum 
regime—, even when interactions are very strong. We evaluate these bounds for harmonic and 
strongly anharmonic junction models and compare with numerical approaches. 
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Does quantum mechanics place bounds on the rate of 
equilibration in a many-body system? When the ini¬ 
tial state of the system describes an inhomogeneity in 
some quantity (e.g., temperature, chemical potential, ve¬ 
locity) , this question can be cast in terms of bounds on 
currents and, in the linear response regime, their corre¬ 
sponding transport coefficients (thermal and charge con¬ 
ductivities, viscosity). Heuristic bounds on such coeffi¬ 
cients based on Heisenberg’s uncertainty principle have 
been conjectured in a wide array of systems, ranging 
from the quark-gluon plasma [1, 2] and ultra-cold atomic 
gases [3, 4], to incoherent “bad metals” [5, 6]. The consid¬ 
erable utility of such bounds is that they can give physi¬ 
cal insight into the nature of transport in regimes where 
interactions are strong and perturbative calculations of 
transport break down [2-7]. Unfortunately, as far as we 
know, none of these bounds have yet been made rigor¬ 
ous [8]. 

Bounds on heat transport are of special value, connect¬ 
ing to foundational questions of information and entropy 
flow [9, 10]. At the same time, there is a pressing need to 
understand the limits of nanoscale devices such as semi¬ 
conductor nanowires [11], carbon nanotubes [12], silicon 
membranes [13], and molecular chains [14-16], to manip¬ 
ulate and transport heat [17-21], an enterprise for which 
interactions may prove crucial [22]. For systems in which 
a heat current arises in a subsystem connecting two ther¬ 
mal baths at temperatures Tl and Tr (see Fig. 1) with 
the current described by the Landauer-type expression 

Jq = ^ J dajhwT{uJ,TL,Tii)[nL{uj) - riRiuj)], ( 1 ) 

it is well known that the heat current is bounded by [9, 
10, 23, 24]: 

jQ<^(r|-r|). (2) 

riR and ur are the Bose (or Fermi) thermal distribu¬ 
tion functions of the left and right thermal baths and 



FIG. 1: Schematic figure of a nanojunction comprising a cen¬ 
tral subsystem described by the hamiltonian Hs, connecting 
two thermal baths at temperatures T^, Tr. A heat current 
Jq arises when Tl ^ Tr. 

T{uj,Tr,Tr) is a generalized transmission function, pos¬ 
sibly including interaction effects. For small temper¬ 
ature differences, (2) reduces to Jq < kq{Tr — Tr), 
where kq = 7TkgT/6h is the quantum of thermal con¬ 
ductance [25]. 

Equation (2) is a weak quantum bound, reflecting only 
unitarity, as can be seen by inserting T(a; > 0, Tr, Tr) < 
1 (for bosons [26]) into (1). This means that (2) will 
only be saturated by systems exhibiting ballistic trans¬ 
port [27-29], for which T(w > 0) = 1 [25]. In nanoscale 
systems showing deviations from ballistic transport, as 
happens when the thermal energy ksiTp + Tr)/2 is 
greater than the energy scale(s) huj* of the low-energy 
degrees of freedom of the subsystem hamiltonian (but 
much smaller than the energy scale hujc at which the 
subsystem behaves classically) or when interactions are 
strong, (2) will not provide a close approximation to the 
actual heat current. 


We emphasize that (1) is not limited to weak in¬ 
teractions in the subsystem separating the two baths. 
As shown by Meir and Wingreen for charge conduc¬ 
tion [30], and subsequently extended by others to heat 
transport [31-35], (1) is an exact expression for the cur¬ 
rent whenever the coupling at the left “lead” is propor¬ 
tional to that at the right [c.f. (5)]. The identification of 
T{uj,Tl,Tr) with nonequilibrium [Tr ^ Tr) correlation 
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functions for the subsystem operators [30-35] allows us 
to relate the heat current to exact spectral sum rules for 
these correlators. 

In this Letter, we derive rigorous quantum mechani¬ 
cal bounds on heat currents described by (1) that pro¬ 
vide a much stronger bound than (2) when hio*/ks ^ 
TltTr <g; Hujc/kB or interactions are strong [this ex¬ 
cludes the ideal gas for which (2) is saturated since it 
is devoid of any energy scale]. Unlike (2), which fur¬ 
nishes a bound in terms of only temperature and funda¬ 
mental constants, our bounds involve a thermodynamic 
quantity-the /-sum rule—that is system-dependent and 
must be evaluated numerically. As with e.g., a conjec¬ 
tured bound on the shear viscosity of quantum fluids 
involving the (system-dependent) entropy [2, 36], it is 
precisely this feature that makes the bound universal, 
since it is saturated by all systems with finite character¬ 
istic low-energy scales at high temperatures. In contrast, 
(2) is only saturated for an ideal gas. Our bounds give 
physical insight into the behaviour of heat transport in a 
wide array of systems while replacing the difficult task of 
calculating time-dependent correlation functions needed 
for T(w, Tl, Tfl) with the much easier task of calculating 
a simple operator expectation value (the sum rule). 

Heat current through a nanojunction coupled linearly 
to two boson heat baths —We focus here on boson (e.g., 
phonon and photon) heat transport through the system 
shown schematically in Fig. 1. The straightforward ex¬ 
tension to electron heat transfer will be considered in 
another publication [37]. We model this system by the 
hamiltonian 


^ = E E 

I'—L.R k 


+ '^9uk^nibik + ^’'k) 

n 




( 3 ) 

Here b^k annihilates a bath particle of energy Huj^k 
with quantum number (usually a momentum) k in the 
v = L/R bath. Their occupancies obey the equilib¬ 
rium Bose distributions n,^{ui^k) = [exp(/ia;i/fc/3iy) — 1]“^, 
= {kBT,y)~^. The subsystem connecting the two 
baths is described by the hamiltonian where 

{Ayi} denotes the set of operators describing the quan¬ 
tum states of the subsystem; ft is in general a composite 
index denoting e.g., the quantum levels for a harmonic os¬ 
cillator junction, the spin index in a spin-boson model, or 
pairs of levels in a more general anharmonic junction (in 
which case one could choose e.g. Aft = lf)(/|, n = {*,/} 
describing the transition between two levels), g^f. is the 
coupling between Afi and the bath degrees of freedom. It 
is characterized by the generalized spectral function [38] 
(not to be confused with the heat current Jq to which 
we always attach the subscript Q) 


j"’" (w) = ^ E - ^.k). (4) 

k 


As shown by several authors using Keldysh meth¬ 
ods [32-35], analogous to the famous Meir-Wingreen ex¬ 
pression for the charge current through a junction [30], 
one can derive an exact expression for the heat current 
corresponding to the model (3), irrespective of the de¬ 
tails of Hs- When the spectral functions on the left and 
right sides are proportional, J"’" (w) = 20 ^ (w), the 

steady-state heat current between the left and right baths 
has the particularly simple Landauer-type form [c.f. (1)] 


Jq = - 


h? 


ay 


_ poo 

/ da;a;Imxs.n'(w)Js,s'(w)[nL(w)-n_R(a;)]. 

do 


( 5 ) 

an are the dimensionless bath-junction couplings, possi¬ 
bly different at the left and right baths, a = az, -I- an, 
and 7 = 4az,a/{/a^. Xn,nA^) = f dsexp(iu}s)xfi,n'(s), 
s = t — t', is the Fourier transform of the retarded corre¬ 
lation function for the operators Af^: 


Xn,n'{.t-t') ^ '-Q{t-t'){[Anit),AlAt')]). ( 6 ) 

As we show in the Supplemental Materials, in steady- 
state, where the density matrix p[t) = p is independent 
of time, Xn,n' obeys the well-known “/-sum rule” [42] 

1 1 - 
= ~ dujujlinxn,n'{<^) = W’An]]), 

( 7 ) 

even for the nonequilibrium situation where the expec¬ 
tation value (• ■ •) = Tr[p(- • •)] in (6) and (7) is not the 
thermal equilibrium one. 

Quantum bound on the heat current —We now make 
use of the /-sum rule to derive a rigorous bound on the 
heat current (5). For simplicity, we will make the (nat¬ 
ural) assumption that the form of the dissipation is the 
same for all of subsystem states, 

where is a dimensionless number that acts as a 

“flag” to enumerate the appropriate form for the op- 
erator(s) Aa^An, coupling the thermal baths. We also 
restrict ourselves to the usual situation where the dissi¬ 
pation is “ohmic” with an exponential cutoff [38]: 

= hfi,fi'^exp{-U!/uJc). ( 8 ) 


In the Supplemental Materials we show how to generalize 
these bounds for other forms of dissipation. 

Using (8) as well as the inequality 0 < hjuj[nL{oj) — 
nfl;(a;)]/A:B(Tz - Tr) < 1 for w > 0 and Tr > 7A, the 
positivity of the change in entropy, requiring Imy^j (w > 
0) > 0 Vw [43], immediately leads to the bounds 

Jq < -^^fiksiTL — Tr) (9) 

n,n' 

and hence, 

— T ^ ^n,n' ■ 

n,fi' 

( 10 ) 
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Equations (9) and (10) are the central results of 
this Letter, giving rigorous bounds on the heat cur¬ 
rent through a nanojunction. Taking 
they admit a simple physical interpretation: 2h^mA.A = 
([^'1', [H, ^]]) is the minimum energy Ej^in that can be de¬ 
posited in the system by the operator A [44] . Our bounds 
can thus be understood as expressing the Heisenberg- 
limited time r > h/Emin for quasiparticles to transport 
heat in Jq ^ ksiTL — TjA/t. This provides a very natu¬ 
ral extension of the quantum of conductance bound which 
instead effectively assumes Elmin ~ + Tr). The 

assumption of Ohmic dissipation is not crucial and any 
form of dissipation would result in sum-rule bounds of 
the form given above (see Supplementary Materials). In 
contrast, the assumption that the left and right spectral 
functions have the same form seems to be crucial: with¬ 
out this, the Meir-Wingreen expression for the current no 
longer involves an integral over a sign-definite function, 
and there is no bound in general. 

For completeness, we give the general (although not 
very transparent) result for the multi-level model 


Hs = '^E^\n){n\; Aij = \i){j\. 

n 

For this model, the /-sum rule is 


yp I Akig{,p 

- ^ ^ )}, 


yp 


( 11 ) 


( 12 ) 


where Aij = {Aij). As in the spin boson and harmonic 
junction models we consider below, it is usually the case 
that the coupling assumes a simple form and one can 
always find a basis where the relevant sum rule is greatly 
simplified, independent of the bath degrees of freedom. 
First, we briefly discuss the regimes where our sum rule 
bounds are saturated. 

As mentioned earlier, a useful feature of our bounds is 
that they are expected to be saturated at high temper¬ 
atures. Specifically, one needs to be at sufficiently large 
temperatures (but smaller than Uc) that 71 ^( 10 )—nfi{u}) ~ 
ksiTL — Tfi)/}Ej over the range of frequencies where 
Imxfj ^j/(a;) is appreciable. This is the case as long as 
both temperatures Tl,T[{ ^ fvjj*/ks are larger than 
the scale(s) huj* that characterizes the central subsys¬ 
tem. For a single-level harmonic junction, w* = ojq is just 
the harmonic oscillator frequency, while in the nonequi¬ 
librium spin-boson model, uj* A^, the renormalized 
tunnelling frequency [38] (related, for ohmic dissipation, 
to the Kondo temperature Tr- [35]). 

Harmonic junction — We now evaluate the sum-rule 
bound for the case of a single-level harmonic junction 

Hs = huJo{a^a + 1/2)] A = a-|-a^. (13) 



FIG. 2: Heat current through a harmonic junction at low (a.) 
and high (b.) temperatures as a function of the temperature 
difference Tl — Tr between the left and right baths for oz, = 
aR — 0.4 (a = 0.8, 7 = 1 ). The dotted line is the exact result, 
the dashed blue line is the ballistic bound (2), while the solid 
red line is the sum-rule bound (15). 


This model is exactly solvable and so, sum rule bounds 
on the current are not needed. However, it is known to 
saturate the “ballistic” bound (2) at low temperatures, 
and so provides a useful system to compare this result 
with our new bound (9). 

Reflecting the fact that this model admits an exact 
solution [33, 45-48], 


/ N _ _ 2u)o/h _ 

- Wq - i7ra;o[JL(a;)-I-Jfl(w)]/2’ 


(14) 


the /-sum rule (7) assumes a simple form, equal to a 
constant: mo,o = wo//i. Therefore, 


jQ<^UJokB{TL-TR). (15) 


Using (14) in (5) with ohmic dissipation, the heat cur¬ 
rent is plotted in Fig. 2 and compared with the ballistic 
bound (2) as well as (15). At low temperatures, the exact 
heat current nearly saturates the ballistic bound. In con¬ 
trast, at higher temperatures, Tr, Tr > fiwo/fcs, the heat 
current is far below the ballistic bound but very nearly 
saturates the sum-rule bound (15). This trend is more 
clearly seen in the thermal conductance, which we plot 
in Fig. 3 together with the thermal conductance bound 
(10) and the quantum of thermal conductance. 


These results confirm the picture laid out earlier: at 
temperatures greater than the characteristic energy scale 
hujf} of the subsystem, the heat current is bounded by coq 
and not the temperature. Hence, the sum-rule bounds 
(9), (10) do far better at high temperatures than the bal¬ 
listic bound (2). Interestingly, because the operator A in 
(13) that couples to the baths is proportional to the dis¬ 
placement operator x of the harmonic oscillator, adding 
anharmonic terms such as V (x) oc x^ to the susbsystem 
hamiltonian will not alter (15). For a minimal model of 
an anharmonic oscillator that gives rise to a nontrivial 
bound, we need to turn to the spin-boson model. 
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FIG. 3: Temperature dependence of the exact thermal con¬ 
ductance K, (dotted line), the quantum of thermal conduc¬ 
tance KQ (blue dashed line), and the quantum sum-rule bound 
7 ra 7 tJofcs /4 (red solid line) for a = 0.8, 7 = 1. At low temper¬ 
atures, the thermal conductance saturates kq, while at high 
temperatures it approaches the sum-rule bound. 


Non-equilibrium spin boson model —We now apply our 
bounds to the spin-boson model, which does not admit 
an exact solution (for a review of various approximation 
methods, see Refs. [49-51]). For this model, 

Hs = ^^x; A = a^, (16) 

hz,z = 1, the /-sum rule is mz,z = and 

hence, 

Jq < ^kBin - Tn) {-A{ax)) (17) 

and 

( 18 ) 

The above results are most easily obtained from niz^z = 
{[az,[H,<Tz]])/ 2 h‘^; they also result from the more gen¬ 
eral (11) and (12) by setting gl^p = g^p[Sip 6 j ,2 + ^i.2<5j,i] 
(corresponding to coupling between the baths and the ax 
operator), HA = E 2 — Ei, and rux^x = (m-12,12 + ^12,21 + 
^21,12+W2i,2i) (which follows from ax = |1)(2|-|-12)(1|), 
and afterwards rotating the spin basis: az ^ ax, ax ^ 
—Oz- The dependence on the bath degrees of freedom 
that arises in (12) cancels out in mz,z- 

In contrast to the exactly solvable harmonic junction 
model, the right-hand sides of these bounds must be 
evaluated numerically. This is straightforward to do ex¬ 
actly using the quasi-adiabatic path-integral (QUAPI) 
approach [52, 53] and in Fig. 4, we plot the sum-rule 
conductance bound (18) for a = 0.2,7 = 1 together 
with the quantum of thermal conductance kq . For com¬ 
parison we also show the thermal conductance taken 
from a recent quantum Monte Carlo (QMC) calcula¬ 
tion [35] as well as the weak-coupling Bloch-Redfield 
(BR) result which, in the limit Wc A is given by [49] 
K = (/i^7r/8fcB)(A^/T^)a7/ sinh(/3fiA). 



FIG. 4: Temperature dependence of the thermal conductance 
bounds (red solid and blue dashed lines) for a = 0 . 2,7 = 1 . 
The QMG value for k [35] (dotted line) saturates the sum- 
rule bound (red solid line) for temperatures greater than the 
Kondo temperature Tk, but is always much less than the 
quantum of thermal conductance (blue dashed line). The 
Bloch-Redfield value (green dot-dashed line) violates the 
bound at intermediate temperatures. The sum-rule bound 
(18) is calculated using QUAPI for ujc = lOA and the numer¬ 
ical algorithm developed in Ref. [53]. 


As expected, the sum-rule bound is saturated from 
below by the (in principle) exact QMC result at tem¬ 
peratures above the Kondo temperature [35, 38] Tk = 
ga{hA/kB){A/u:c)°'/^^~°'1 for a < 1 and 0 for a > 1. 
ga = [r(l — 2a) cos(7ra)]^/^^^““\ with F the Gamma 
function. For T ^ Tk, both k [35] and {ax) [54] scale 
as {T/Tk)^°‘~^, as required by our sum rule bound. The 
BR thermal conductance violates the sum rule bound at 
intermediate temperatures T ~ Tk- For higher tempera¬ 
tures, it scales as 1/T and lies below the sum-rule bound 
as well as the QMC result. At weaker coupling, a < 0.1, 
the BR value is below the bound for all temperatures. 

The sum rule bounds (17) and (18) show the crucial 
role played by the off-diagonal coherence (ax) in the heat 
current. For small coupling a, (ax) is independent of a 
and the heat current will increase linearly with the cou¬ 
pling. Ultimately, however, for values of this coupling 
beyond the so-called Toulouse point a = 0.5 [38, 55], the 
spin subsystem is strongly entangled with the environ¬ 
ment, the coherence (ax) is greatly suppressed, and the 
thermal conductance becomes asymptotically small [35]. 


Conclusions —We have derived rigorous quantum 
bounds on heat transport through nanojunctions that 
extend the well-known quantum of thermal conductance 
bound to better describe interacting systems. These 
bounds are universal in the sense that for any system 
with characteristic low-energy energy scales, the actual 
heat current saturates our bounds for temperatures above 
these scales. In this way, our bounds will always outper¬ 
form the ballistic bound at these temperatures, while at 
temperatures below these scales, the ballistic bound is 
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generally stronger [see Figs. 3 and 4]. Ideal gases that 
saturate the ballistic bound are devoid of such scales 
and our bound will do worse (it is formally divergent!). 
For all other systems with energy scales introduced by 
e.g. interactions or an ultraviolet cutoff (as in few-level 
systems), our bounds are well-saturated at intermediate 
temperatures, opening the door to investigating room- 
temperature heat transport in nanoscale devices when 
interactions are strong. It will be particularly interesting 
to revisit foundational questions of entropy flow [9, 10] 
when interactions are strong and also to study how in¬ 
teractions can be used to engineer quantum heat devices 
such as thermal diodes [56]. 
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SUPPLEMENTAL MATERIAL 


Spectral properties of the nonequilibrium correlation 
function — Here we discuss the spectral properties of the 
in-general nonequilibrium retarded correlation function 
for the operators A and B [note that for brevity, we 
adopt a slightly different subscript notation in Eq. (6) 
in the main text]: 

XAMt, t') ^ - t')([i(t), B\t')]). (19) 

The expectation value that enters this expression can be 
expressed in terms of the eigenstates | a), 16) of the hamil- 
tonian H as 


{Q{t)) = ^Pab{a\Q{t)\b) 

ab 




( 20 ) 


Here e.g., Ea is the eigenvalue corresponding to the 
eigenstate |a): H\a) = Ea\a). pab are the matrix ele¬ 
ments of the in general time-dependent density matrix 
by P = J2ab Pa-b\b){a\; in thermal equilibrium, it as¬ 
sumes the usual form pab = Sa,bexp{—PEa)/Z, where 
Z = Tr[exp(—/3i7)] is the partition function. We will 
generically be interested in the nonequilibrium situation, 
where Tl A Tr. As with the equilibrium situation, how¬ 
ever, in steady-state, where the left-hand-side of (20) is 
time-independent, the summation in this expression col¬ 
lapses and one can replace pab with PaaSa,b- 

Using (20), one can write down the spectral representa¬ 
tion for xa,b using the exact eigenstates and eigenvalues 
of the full hamiltonian. Doing this, one obtains 


Xa,b{^A') = 


eh'' 


-Q{t-t')^Pac 

abc 




{a\B^b}{b\A\c) 


( 21 ) 


We suppose that the system settles down to a steady- 
state configuration after some characteristic time t*. For 
t, t' both much larger than this time scale, the two-time 
retarded correlator XAA{t, t') only depends on the differ¬ 
ence between t and tb 


XA.Bit,t') ^ XA,B{t-t'). (22) 


This means that the summation over the indices a, &, c 
in (21) collapses, with c = a. Fourier transforming the 
resulting expression gives 


Xa,b{uj) = 

1 ^ r (a|Ht|6)(6|i|a) 

hLO + Eba+ i0+ 

ab L 


(a|i|b)(b|Ht|a) ' 
hw- Eba+ *0+ ’ 


(23) 


where E^a = Eh — Ea- The imaginary part that enters 
the expression for the heat current is 

ImxA,B(w) = ^'^{paa-Pbb){a\A\b){b\B^a)S{fiJjj-Eba). 

(24) 

Multiplying (24) by w and integrating yields 


roo 2 _^ 

- / dwa;ImXA,B(w) = ■^'^{Paa - Pbb)EbaX 

^ ” ab 

{a\A\b){b\B^a) = ( 25 ) 


Making use of the fact that ImxA,B(—w) = “IniXA.B(w), 
setting A = Afi,B = and recalling that pab can be 
replaced by PaaSa,b in ( 20 ) in steady-state gives the /- 
sum rule shown by Eq. (7) in the main text. 

Other forms of dissipation —For simplicity, we assumed 
ohmic dissipation in the main text. Here we show how 
to extend our results to other forms of dissipation. 

Any form of dissipation Jfi,n'i^) leads to a 
bound on the heat current: one only needs to 
find a function A Jnfi'iyi) such that 

Jjj dujEfi^fiAu!)lmXfi,n'i^) can be calculated (this can al¬ 
ways be done). For F{uj) = uj\ where I is an inte¬ 
ger and we have dropped the n,n' subscript for clarity, 
such sum rules are straightforwardly given as a commu¬ 
tator expression, as in (25). In all cases of current in¬ 
terest, including “super-ohmic” dissipation with / > 1 , 
Lorentzian, and Gaussian forms, the spectral function 
can be bounded by such a simple polynomial expression. 
At the same time, for e.g., a Lorentzian spectral func¬ 
tions of the form J(w) = a;^[l -I- 77 ( 0 ; — ujr)‘^]~^ < 
improved bounds can be obtained by directly evaluat¬ 
ing dujJ{uj)linXfi,n'{^)- This can be done by apply¬ 
ing the operator-product expansion of Kadanoff, Wilson, 
and Polyakov [1]. (See Refs. [2, 3] for examples of such 
“smeared” sum rule calculations.) The same technique 
can be used to calculate sum rules for non-analytic E{uj), 
including power-law forms with non-integer exponents. 

For these more general forms of dissipation, in addi¬ 
tion to the high-temperature criterion outlined in the 
main text, in order for the exact heat current to satu¬ 
rate the resulting sum rule bound, one needs the spec¬ 
tral function to be approximately power-law over the 
range of frequencies where Imxn.n' is appreciable. For 
the Lorentzian form discussed above, this means that 
Imxn,n' must be peaked in the vicinity of the resonance 
frequency ujr. Even if this condition is not satisfied, 
we emphasize that a.) the bound remains rigorous and 
b.) one can derive improved sum rule bounds along the 
lines suggested above, and so at sufficiently large tem¬ 
peratures, hio* <C ksTL^ksTR HlOc, the appropriate 
sum-rule bound will always be very nearly saturated. 
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